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Abstract 

Consider the Quadratic Assignment Problem (QAP): given two matrices A and 
D, minimize { trace (AXDX T ) : X is a permutation matrix }. New lower bounds 
were obtained recently (Mittelmann and Peng [5]) for the QAP where D is either 
the Manhattan distance matrix of a rectangular grid or the Hamming distance 
matrix of a hypercube. In this note, we show that the results in [5J lllj extend 
to the case where D is a spherical Euclidean distance matrix, which includes the 
Manhattan distance matrix and the Hamming distance matrix as special cases. 

1 Introduction 

The Quadratic Assignment Problem (QAP)( see, e.g., [4] and references therein) is the 
problem: 

min tr&ce(AXDX T ), 
xen 

where A and D are two given matrices of order n and n is the set of n x n permutation 
matrices. New lower bounds were obtained recently (Mittelmann and Peng [Sj) for the 
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QAP where D is either the Manhattan distance matrix of a rectangular grid or the 
Hamming distance matrix of a hypercube. Mittelmann and Peng exploit the fact that 
for such matrices D there exists a scalar A such that XE — D is positive semidefinite, 
where E is the matrix of all l's. In this note, we show that the results in [Hill] extend 
to the case where D is a spherical Euclidean distance matrix (EDM), which includes 
the Manhattan distance matrix and the Hamming distance matrix as special cases. 

1.1 Notation 

The positive semidefiniteness of a symmetric real matrix A is denoted by A y 0. The 
set of symmetric real matrices of order n is denoted by <S n . J n denotes the identity 
matrix of order n and e n denotes the vector of all l's in W 1 . E n denotes the n x n 
matrix of all l's, i.e., E n = e n e^. The subscript n will be omitted if the dimension of 
J n , e n and E n is clear from the context. A^ denotes the Moore- Penrose inverse of A. 
For two matrices A and B, A®B denotes the Kronecker product of A and B. For a real 
number x, \x\ denote the ceiling of x. Finally, diag(A) denotes the vector consisting 
of the diagonal entries of a matrix A. 

2 Preliminaries 

The following well-known [5] lemmas will be needed in some of the proofs below. 

Lemma 2.1 Let A and B be two n x n real matrices. 

1. (A T y= (AY := AV. 

2. (AA T )^=A^A^. 

3. If A T B = and BA T = 0, then (A + B)t =A*+B*. 

4. If A T B = and BA T = 0, then rank (A + B) = rank A + rank B. 

5. If A has full column rank then A^ = (A T A)~ l A T . 

6. Let C and Q be two real matrices of orders k x k and n x k respectively, where 
Q T Q = I k . Then (QCQ T ) t = QC^Q T . 

Lemma 2.2 Let A and B be two n x n real matrices. 

1. (A = At (g>5t. 

2. Let Ah and Bt.0. Then (A®B)t_ 0. 
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3. rank(Aig)B) = rank (A) rank (B). 

Lemma 2.3 (Generalized Schur Complement) Given the real symmetric parti- 
tioned matrix 

A B ' 
B T C 



M 



Then M y if and only if C y 0, A — BC^B T y 0, and the null space of C is a subset 
of the null space of B . 

3 Euclidean Distance Matrices (EDMs) 

An n x n matrix D = (dij) is called a Euclidean distance matrix (EDM) if there exist 
points p 1 , . . . , p n in some Euclidean space W such that 

dij = \\p l ~ p'f for all i, j = 1, . . . , n, 

where || || denotes the Euclidean norm. Moreover, the dimension of the affine span of 
the points p 1 , . . . ,p n is called the embedding dimension of D. Without loss of generality, 
we make the following assumption. 

Assumption 3.1 The origin coincides with the centroid of the points p , . . . ,p n . 



It is well known |12t uA\ that a symmetric matrix D whose diagonal entries are all O's 
is an EDM if and only if D is negative semidefinite on the subspace 

M := {x G R n : e T x = 0}, 

where e is the vector of all l's in R n . It easily follows that the orthogonal projection 
on M is given by 

J -I- ee T /n = 1- E/n. (1) 

Let S n denote the set of symmetric n x n real matrices and let Sh and Sc be two 
subspaces of S n such that: 

Sh = {A G S n : diag(A) = 0} and S c = {A E S n : Ae = 0}. 

Furthermore, let [5] T : Sh — > Sc and K. : Sc — > Sh be the two linear maps defined by 

T(D) := ~JDJ, (2) 
K{B) := diag(S) e T + e {di&g{B)) T -2B. (3) 
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Then it immediately follows that T and K, are mutually inverse between the two 
subspaces Su and Sc\ and that D is an EDM of embedding dimension r if and only if 
the matrix T{D) is positive semidefinite of rank r [5]. Moreover, it is not difficult to 
show (see, e.g., [U EJ [13]) that if D is an EDM of embedding dimension r then rank 
(-D) is equal to either r + 1 or r + 2. 

Note that if D is an EDM of embedding dimension r then T{D) is the Gram matrix 
of the points p 1 , . . . ,p n , i.e., T(D) = PP T , where P is n x r and p l is the ith row of 
P. i.e., 

P = 



lT 

P 



(4) 



P is called a configuration matrix of D. Also, note that P T e = which is consistent 
with our assumption that the origin coincides with the centroid of the points p 1 , . . . ,p n . 

3.1 Spherical EDMs 

An EDM D is called spherical if the points p 1 ,]? 2 , . . . that generate D lie on a 
hypersphere. Otherwise, D is called non- spherical. The following characterization of 
spherical EDMs is well known. 

Theorem 3.1 ( [10|, [6], I13j ) Lei D ^ be a given nx n EDM with embedding dimen- 
sion r, and let P be a configuration matrix of D such that P T e = 0. Then the following 
statements are equivalent. 

1. D is a spherical EDM. 

2. Rank (D) =r + l. 

3. The matrix A ee T — D is positive semidefinite for some scalar A. 
4- There exists a vector a in W such that: 

Pa= l -Jdiag{T(D)), (5) 

where J is as defined in (CP. 

Two remarks are in order here. First, spherical EDMs can also be characterized in 
terms of Gale transform [2] which lies outside the scope of this paper. Second, if D is 
an n x n spherical EDM then the points that generate D lie on a hypersphere of center 
a and radius p = (a T a + e T De/2n 2 ) 1 ^ 2 , where a is given in ([5]) (see [E]). Furthermore, 
as the next theorem shows, the minimum value of A such that XE — D >z is closely 
related to p. 
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Theorem 3.2 ( \10[ I13| ) Let D be a spherical EDM generated by points that lie on a 
hypersphere of radius p. Then X* = 2p 2 is the minimum value of X such that XE — D y 
0. 

The following lemma gives an explicit expression for the radius p which is more conve- 
nient for our purposes. 

Lemma 3.1 Let D be a spherical EDM generated by points that lie on a hypersphere 
of radius p. Then 

2 _ fDe e T D(T(D))^De 
9 2n 2 + 4n 2 ' 1 ' 

Proof. Let D be a spherical EDM and let Q = [-^U V] be an orthogonal matrix, 
i.e., J = VV T . Then XE - D y if and only if Q T (XE - D)Q h 0. But 



Q T (XE - D)Q 



An — e De/n —e DV/y/n 
-V T Dej4n -V T DV 



Clearly -V T DV = 2V T T(D)V y 0. Thus it follows from Lemma E3] that XE — D y 
if and only if An - e T De/n - e 1 " D\ r (V T T{D)V^V 7 "Dej2n > and the null space of 
V T DV is a subset of the null space e T DV . But it follows from Property ([6]) of Lemma 
EI that V{V T T{D)V)W T = {VV T T{D)VV T )^ = T{D)^ since VV T = J and since 
(T(D))e = 0. Moreover, since D is a spherical EDM, it follows that the null space 
of V T DV is a subset of the null space e T DV [2]. Thus, XE — D y if and only if 
An — e 1 De/n — e T D(7~(D))^ De/2n > and the result follows by Theorem 13.21 

□ 

Note that (T(L>)) t =(PP T ) t = P tT P f . Thus it follows from © that p 2 = a T a + 
e T De/2n 2 = e T De/2n 2 + e T DP^ T Pt De/An 2 . Therefore, 

a=l^De = l(P-Pr^De, (7) 

where P is a configuration matrix of D. 

An interesting class of spherical EDMs is that of regular EDMs. A spherical EDM 
D is said to be regular (also called EDM of strength 1 [9l [TO]) if D is generated by 
points that lie on a hypersphere centered at the centroid of these points [7]. i.e., 

na = P T e = 0, 

since we assume that the centroid of the points p ,...,p n is located at the origin. 
Therefore, if D is a regular EDM then diag(7~(D)) = p 2 e. Hence, it follows from ([3]) 
that 

De = 2np 2 e, 

where p 2 = e T De/2n 2 . Thus we have the following characterization of regular EDMs 
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Theorem 3.3 ([.3 13 US]) Let D be an EDM. Then D is regular if and only if e is 
an eigenvector of D. 

4 Main Results 

Theorem 4.1 Let D\ be an m x m EDM of embedding dimension r\, and let D2 be an 
n x n EDM of embedding dimension r 2 ■ Then D = E m (g> D2 + D\ (g) E n is an EDM of 
embedding dimension r\ + r 2; where E m is the m x m matrix of all 1's and <S> denotes 
the Kronecker product. 

Proof. I nm = I m ®I n and E nm = E m ® E n . Thus 

T(E m ®D 2 ) = ~(I m ®I n -—E m ®E n )(E m ®D 2 )(I m ®I n -—E m ®E n ) 
2 nm nm 

111 

= —-Em ® {In E n )D2(I n E n ) 

Inn 
= E m ® T(D 2 ) h 0. 

Similarly T(D 1 ® E n ) = T{D X ) ®E n yO. Therefore, T{D) = E m ® T{D 2 ) + T(D X ) ® 
E n y 0. Hence D is EDM. 

On the other hand, we have from Lemma 12.11 that rank (T(D)) = rank (E m (8) 
T(D 2 ) + T(-Di) E n ) = rank (E m T{D 2 )) + rank (T(L>i) ® E n ) = rank (T(£> 2 )) + 
rank (T(-Di)) = r 2 + n. Thus the embedding dimension of D is equal to r\ + r 2 . 

□ 

Theorem 4.2 Let D\ be a spherical EDM of order m generated by points that lie on a 
hypersphere of radius p\ , and let D2 be a spherical EDM of order n generated by points 
that lie on a hypersphere of radius p2- Then D = E m <g> D2 + D\ ® E n is a spherical 
EDM generated by points that lie on a hypersphere of radius p = {p\ + p^) 1 ^ 2 ■ 

Proof. XE nm — D = E m ® (\ 2 E n — D 2 ) + (XiE m — D x ) ®E n where A = X x + A 2 . 
Since D\ and D 2 are spherical EDMs it follows that A 2 -B n — D2 y for A 2 = 2p 2 , and 
\ x E n -D x yQ for Ai = 2p\ . Therefore, for \ = 2p\ + 2p\ we have XE nm — D y and 
hence D is a spherical EDM generated by points that lie on a hypersphere of radius 
P < (Pi + pl) l/2 - Next we show that p = (pj + pl) 1/2 . 
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By Lemma 13. II we have 

„2 



e T De e T D(T(D)^De 
2{nm) 2 4(nm) 2 

16m 



eT,D2e n , eLDie 



2n 2 2m 2 

+{(mei (8) e£D 2 + e^I>i ne T n )(E m ® T(D 2 ) + T(£>i) ® £n) f - 

...(men, (8) Z?2e n + -Die m ® ne n )}/4(nm) 2 
e^D 2 e n e^Die m 



2n 2 2m 2 

+IK ® e^ 2 + e&Z* ® nel)(^ ® (T(L> 2 ))t + (T(Di))t ® % 
...{me m ® -D 2 e n + £>ie m <g> ne„)}/4(nm) 2 

2re 2 + 2m 2 + { e n D 2 (T(I>2)) T £>2e„ 
+n 2 e^D 1 (T(Z?i)) t I?ie m }/(4(nm) 2 
P2+P1 



□ 

Example 4.1 Xei G n = (gy) 6e i/ie Manhattan distance matrix of a rectangular grid 
consisting of 1 row and n columns. Then 

9ij = 

Let q 1 , . . . ,q n be the points in M n_1 such that the first i — 1 coordinates of q 1 are l's 
and the remaining n — i coordinates are O's. Then it immediately follows that q , . . . , q n 
generate G n and q 1 , . . . ,q n lie on a hypersphere of radius p = \{n — 1) 1//2 and centered 
atb = (1/2, 1/2,... ,1/2) T . ThusG n is a spherical EDM of embedding dimension n — 1. 
Note that in this example, we don't assume that the origin coincides with the centroid 
of the points q , . . . , q n in order to keep the expressions of the q % 's and b simple. 

Now consider a rectangular grid of m rows and n columns. Let jy be the Manhat- 
tan distance between the grid point of coordinates (i, j) and the grid point of coordinates 
(k,l). Then 

dij,ki = \i-k\ + \ j - l\. 

Let s = i + n{j — 1) for j = 1, . . . , m and i = 1, . . . ,n; and let t = k + n(l — 1) for 
I = 1, ... ,m and k = 1, . . . , n. Then 

j = \s/n\ , i = s — n{ \s/n\ — 1) and I = \t/n\,k = t — n{\t/n\ — 1). 
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Define the nm x nm matrix D = (d s t) such that d s t = ckj ki- Then it follows [8] 
that 

D = E rn ® G n + G m <g> E n , (8) 

where G n is as defined in Example 14.11 Equation (JBJ) follows since (E m ®G n ) s t = \i— k\ 
where i = s — n(\s/n\ — 1) and k = t — n{\t/n\ — 1); and since (G m ® E n )ki = \j — l\ 
where j = \s/n\ and I = \t/n\. 

Thus we have the following two corollaries. 

Corollary 4.1 Let D be the mn x mn Manhattan distance matrix of a rectangular 
grid of m rows and n columns. Then D is a spherical EDM of embedding dimension 
n + m — 2. Furthermore, the points that generate D lie on a hypersphere of radius 
p = \{n + m-2) 1 / 2 . 

Proof. This follows from Theorems 14.11 and 14.21 

□ 

Corollary 4.2 (Mittlemann and Peng [8, Theorem 2.6]) Let D be themnxmn 
Manhattan distance matrix of a rectangular grid of m rows and n columns. Then 
(n + m- 2)E mn /2 - D >z 0. 

Example 4.2 Consider the Hamming distance matrix of the r-dimensional hypercube 
Q r whose vertices are the 2 r points in W with coordinates equal to 1 or 0. Letp 1 , . . . ,p 2 ' 
be the vertices of Q r and let D = (rfy) be the 2 r x 2 r matrix where dij is the Ham- 
ming distance between p l and pi ; i.e., dij = Y?k=l \Pk ~ Pk\' Then dij is also equal 
t° Sfc=i(^A: ~ Pk) 2 = \P l ~P*\\ 2 - Therefore, D is an EDM. Furthermore, the points 
p 1 , ■ ■ ■ ,p 2T lie on a hypersphere, centered at their centoid, of radius p = \r 1 / 2 . Thus D 
is a regular EDM of embedding dimension r. 
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